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\Q ' Abstract 



The system of oscillator interacting with vacuum is considered as a problem of 
random motion of quantum reactive harmonic oscillator (QRHO). It is formulated in 
terms of a wave functional regarded as complex probability process ^f s t c (x,t | W (t)) 

£f~) . in the extended space E = i? 1 CEDi?{vr(i)}- This wave functional obeys some stochastic 

differential equation (SDE). Based on the nonlinear Langevin type SDE of second 
order, introduced in the functional space Riw(t)}> the variables in original equation 

f-- are separated. The general measure in the space R^w^y of the Fokker-Plank type 

is obtained and expression for total wave function (wave mixture) ^ br (n;x,t) of 
random QRHO is constructed as functional expansion over the stochastic basis set 
\I/^ C (n; x, t | W(t)). The pertinent transition matrix S br is constructed. For Wiener 
type processes W{t) the exact representation for " vacuum- vacuum" transition prob- 
ability AqI^q is obtained. The thermodynamics of vacuum is described in detail for 
fti the asymptotic space E s = R 1 (£> R{ Ws y. The exact values for Energy, shift and 

expansion of ground state of oscillator and its Entropy are calculated. 
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H ' 1. Introduction 
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All the processes, described by standard quantum mechanical approach, are stochastic 
processes from the point of view of classical dynamics. The natural equivalence between 
Schrodinger and Fokker-Plank equations was used for formulation of quantum mechanics 
as stochastic theory [1], and the procedure of quantization was introduced [2], that takes 
into account the influence of stochastic processes on dynamics. For solution of quantum 
problems different numerical algorithms were proposed for stochastic dynamics (see [3]). 
Note, that in all above approaches the formulation of the main quantum object, that 
is the wave function, was deterministic. We must underline, that deterministic features 
of the physical theory are the outcome of the symmetry of its main equations with the 
change of the sign of time evolution. 

At same time there is a lot of evidences for quantum deterministic description violation 
both in physics (see [4] ) and in chemistry [5-6]. 

In several papers of the authors [7-9], nonstationary multichannel scattering in collinear 
three-body system was formulated as a problem of wave packet evolution in a system of 



body fixed reper, that makes in general case complex, some times chaotic, motion on the 
induced Riemann manyfold. It was shown, that for three-body system there exist an " in- 
ternal time" describing the evolution of the system, and in which the equation of motion 
is not symmetric with the change of sign of time. It means, that at certain conditions the 
wave function can be the object of probability description. 

In present paper we propose a simple, but nontrivial problem of stochastic quantum 
mechanics - the problem of QRHO under Brownian motion. It was shown by the au- 
thors [9], that such a model can correspond for example to the description of bimolecular 
chemical reaction, that goes via the resonance complex. 

2. Description of the problem 

In the case of the random QRHO the equation for the wave function can be written in 
following form 

i5 t ^stc = H (x, t\W(t)) ^stc, -oo <x,t< +oo, 

(2.1) 
H(x,t\W(t)) = \[-d 2 x + Q 2 (t\W(t))x 2 ], d 2 x = d 2 /dx 2 , h=l. 

with the frequency Q (t\W(t)) , and the wave state fy stc (x, t\W(t)) being the functional 
of, in general, complex Markovian process W(t). We shall denote by St the total derivative 
in view of process ^ stc (x, t\W(t)) (see (3.2) ). We shall suppose also, that the frequency 
and wave functional are subjected to following boundary conditions 

lim n(t\W(t))=Q m{out) >0, (2.2) 

lim y stc (x,t\W(t))= i lim d x V 8tc (x,t\W(t))=0. (2.3) 

\x\— >+oo |a:|— >+oo 

In particular case of frequency, being the regular function of "£" , that is Q (t\W(t)) = 
Q (£), the equation (2.1) with initial condition (2.2) has exact solution (see [9], [10]). 

So in our new approach the equation (2.1) is SDE for complex stochastic process 
\I/J C (n; x, t\W(t)), determined in the extended space S = R 1 <g> R^ W (t)}- 

It should be noted that asymptotic behavior of functional ty~^ tc (n;x,t\W(t)) for the 
time t — »■ — oo according to (2.2) have a following form 



V + (n;x,t\W(t))^ V in (n;x,t) = 

t — > — OQ 

~- : — exp{-i{n+^)Q m T-^Q m x 2 }H n (^/n~x 



(2.4) 



2 n n! 

Our main problems are: 

a) to find the conditions on Markovian process W(t), for which the variables in equa- 
tion (2.1) are separated and so the detailed solution ^/ stc (x,t\W(t)) is found; 

b) the evaluation of evolution of average wave function 



K(n;x,t) = (Vf tc (n;x,t\W(t))) {w(t)} , (2.5) 

that describes the state of random QRHO (where (■■■)i W ( t )\ denote the functional inte- 
gration over the total Fokker-Plank measure including the integration over the distribution 
of stationary process W s = lim W (£)); 



c) computation of corresponding transition S ,6r -matrix, and representation for "vacuum- 
vacuum" transition probability Aq!^ for a case of Wiener type process; 

d) investigation of vacuum thermodynamics for the asymptotic space S s = R 1 (£)R{w a } 
(calculation of shift and width of ground state energy, Entropy and Free energy of the 
oscillator, interacting with vacuum). 

3. Solution of SDE for complex process-wave functional ^ stc (x, t\W(t)) 

Let us start from the equation of classical oscillator under Brownian motion 

Z + n 2 (t\W(t))£ = 0, £ = 8 t S(t\W(t)), (3.1) 

with boundary condition (2.2). Note, that dot over functional £ (t\W(t)) is total derivative 
of Ito type 

£ = Stt (t\W(t)) = %£ + U 2 w i + (SwO d t W(t), (3.2) 

where 5y/ -stands for functional derivative 

8 w Z={8Z(t\W(t))/5W(t)} t=t ,. (3.3) 

Like the following from (2.2) condition, the solution (3.1) have a following asymptotic 
behavior 

£(t\W(t)) ^ exp(in in t) (3.4) 

t — > — oo 

Theorem : If the etalon SDE (3.1) take place, then the SDE (2.1) for complex process 
have an exact solution 

M/+ c (n;x,t|^(t))-r^Z^!l 1/2 



x 
f ' ' ■ 1 (3 ' 5) 



2"ra!|£| 



that due to condition (3.4) goes to asymptotic state (2.4) in the limit t — *■ — oo. 

The Prove. 

The solution of equation (3.1) can be represented in the following form 

£ (t\W(t)) = a (t\W(t)) exp [ir (t\W(t))} , a (t\W(t)) = \£ (t\W(t))\ . (3.6) 

It is obvious, that differentials of Ito type exist also for a(t\W(t)) and r(t\W(t)) 
functionals. 

For further analytic investigation of the problem it is useful to introduce the scales 
of length a (t\W{t)) and time r = r (t\W(t)) /fii n . How one can see, this scales have a 
stochastic character unlike the case of regular problem of parametric quantum oscillator. 

Going to investigation of SDE (2.1) let us make transformation x — > y = xja (t\W(t)) , 
then equation (2.1) will as follows: 

L(y,t\W(t))V 8t c = 0, *rtc(v,t\W(t)) = V 8t c(x,t\W(t)), (3.7) 



L (y, t\W{t)) = iS t - i-y6 y + -L«SJ - ^-Q 2 (t\W(t)) y 2 . (3.8) 

Representing the solution of equation (3.7) in a following form 

**(v,t\w{t))={ PL g( ^ (t) y jyj ) *^/^w*»J (3 - 9) 

and after transformations t ^ t = r {t\W(t)) /Vt in we get from (3.7)-(3.8): 
2 [A - (aa) /2] ($ + 2y^$) + 2 (m 2 /Q in ) <5 T $ = 



-1/2 {^ - a 2 [2A - 4(7(7"^ + 4<j- 2 A 2 + a 2 Q 2 (t\W(t))] y 2 } $ 



(3.10) 



Thus after transformation (x, t) — > (y, r) and substitution (3.9) from (2.1) we arrive to 
equation (3.10), where the functionals a(t\W(t)), r(t\W(t)) and A(t\W(t)) still remain 
to be determined. For their determination let us subject them to realization of following 
conditions: 

r(t\W(t)) = n m /a 2 (t\W(t)), (3.11) 

A (t\W(t)) = & (t\W(t)) a (t\W(t)) /2, (3.12) 

2A - 4(70- _1 A + Aa' 2 A 2 + a 2 Vt 2 (t\W(t)) = Vl 2 Jg 2 . (3.13) 

If we assume that first variation of r (t\W(t)) functional due to W(t) process is equal 
to zero, i.e. Swf (t\W(t)) = 0, then from equation (3.11) it follows, that stochastic time r 
is coupled with natural parameter (usual time) t via the following integral transformation: 

f dt' , 

T = J^WWWY (3 ' 14) 

—00 

As to equations (3.12) and (3.13) it is easy to show, that their combination bring 
the equation (3.1) for complex process C,(t\W(t)). By taking into account expressions 
(3.11)-(3.13) from (3.10) one can obtain the following equation: 

Lo(y,r)$(y,T) = 0, (3.15) 

L (y,r) = i6 T + l -5 2 y - l -n 2 n y 2 . (3.16) 

It is clear that equations (3.15) and (3.16) describe autonomic quantum system, but 
on the stochastic space-time continuum. Solution of (3. 15)- (3. 16) have the following form: 



$(y,r) 



(^A) 1/2 i 1/2 
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exp l-i(n+-)n in r - -Q in y 2 iH n (Jvt in y\ (3.17) 



Combining (3.9) and (3.17) for the complex process fyf tc (n;x,t\W(t)) one can obtain 
final expression (3.5), that had to be proved. It is clear from (3.5), that had stochastic 
process \&^ c (n;x,t\W(t)) at the limit t — > — oo goes to asymptotic state (2.4). 

In conclusion let us pay attention to the following important feature of complex 
stochastic process ^f tc (n;x,t\W(t)): 

oo 

(vf tc (m;x,t\W(t))*f tc (n; x, t\W(t))) x = S mn , (...) x =J ... dx, (3.18) 

— oo 

that shows the fact, that wave functionals make up the full ortho normal basis. 

4. The derivation of Langevin equation for the real stochastic 
process 9{t) 

Now, after determination of the basis in the space of complex functionals \P f tc (n; x, t\W(t)), 
we can pass to construction of expression for averaged wave state ^> br (n; x, t) of quantum 
random oscillator. For this purpose at first it is necessary to determine the measure of 
functional space R{w{t)} on which stochastic process \&J C (n;x, t\W(t)) will be averaged. 
Returning to equation (3.1) let us note, that in general case its analysis is very diffi- 
cult and for its further analytical investigation it is necessary to finalize some features of 

t(t\w(t)). 

Theorem : If the functional £ (t\W (t)) is subjected to the conditions 

SwC (t\W(t)) = 0, 5 W {d t C (t\W(t))} ± 0, (4.1) 

then the stochastic equation (3.1) turns in to nonlinear equation of Langevin type 

e + e 2 + n 2 (t) + F(t\w(t)) = o, (4.2) 

where F (t\W (t)) is the generator of stochastic force. 

The Prove: 

The solution of model equation (3.1) can be represented in the following form: 

t(t\W(t))=So(t)aq>(f *(t'\W(t'))dt'Y (4.3) 

where £o (t) is the solution of equation (3.1) with regular frequency f2 (t). 
After substitution (4.1) into (3.1) and taking into account equation (3.2) one gets for 
$ (t\W(t)) the stochastic nonlinear equation of Langevin type 

$ + 2£ £o _1 $ + $ 2 + F (t\W(t)) = 0, 

(4.4) 
tt 2 (t\W(t))=ttl(t) + F(t\W(t)). 

After transformation 

$ (t\W(t)) = 9 (t\W(t)) - i (t) /Co (t) (4.5) 

we pass from (4.4) to the equation (4.2). Let us note, that transformation (4.5) is 
equivalent to transition to regular moving coordinate system in complex functional space 



i?/$( t )\. As to 9(t\W(t)) functional, it belongs to real functional space Rjqm\. Thus, the 
theorem is proved. 

5. Investigation of Fokker-Plank equation. Determination of the 
measure of functional space R{e(t)} 

Let us pass to derivation of the evolutional equation for condition probability P (9, t\6', £'). 
We shall study the functional of the form 

p(e,t\e',t') = (s[0(t)-9(t'))) {wm , (5.i) 

where 9 (t) is the solution of nonlinear Langevin equation (4.2). After differentiating 
(5.1) over the time and using (4.2) one can obtain 

d t P (9,t\9',t') = -d e (95[9 (t) - «(*')]) = 

(5.2) 

= de{[9 2 + ni(t)\p(e,t\e',t') + (F(t\w(t))S[8(t)-e(t>)]) {wm }. 

The second member in the rhs of equation (5.2) still remains undetermined. For its 
calculation it is necessary to definite the type of stochastic force generator F(t\W(t)). 
As in most interesting cases the F (t\W (£)) = F (i) functional is the gaussian function, 
that in considered problem changes more quickly than £ (t\W (£)), the choice of the model 
of "white noise" for stochastic is quite suitable 



F(t)F(t)) = 2eS(t-t), (F(t)) = 0, e > 0. (5.3) 

Now using the Vick theorem (see [11]) 

(F(t)N(F(t))) {Fm = ^(%) , (5-4) 

\ 0r I {F(t)} 

where N (F (£)) is arbitrary functional of F (t), one can write the following expression 

(FS [9 (t) - 9 (f )]> = -2 ((59 (t) ISF (t)) d e 5 [9 (t) - 9 (t')]) {w(t)} = 

(5.5) 
= -2d e ((59 (t) ISF (t)) 5 [9 (t) - 9 (t')]) {w(t)} . 

Variational derivative of 9 (t) due to stochastic force F (t) equals to e ■ sgn (t — t') + 
O (t — t'). After regularization by standard procedure (in sense of Fourie decomposition) 
one can find it value for t — t'\ e • sgn (0) = \e. Taking into account the above said nota- 
tions now we can obtain now the final expression for Fokker-Plank equation for conditional 
probability: 

d t P (9, t\9', If) = d e { [9 2 + Q 2 (£)] + ed e } P (9, t\0' , t) . (5.6) 

Note, that (5.6) determines the diflusional process, for which 9 (t) is continuous. 



Let the probability be subjected to boundary condition P(6,t\8 , t) = 6(9 — 9), then, 
for small time intervals the solution of equation (5.6) is straightforward [12]: 



P(6,t\6',t') = (27TeAty 1/2 exp\-± ^ ' M, t = t + At. (5.7) 



' 2 ' n 2 (t))At] 2 



2eAt 

It is clear, that the evolution of the system in the functional space R{e(t)} is governed by 
the regular shift with the speed (9 2 + fi^ (£)) modulated by quantum Gaussian fluctuations 
with constant correlations e. 

After this we can establish some properties of the trajectory 9 (t) in the space Rjqm\. 

It is given by the formula (see [12]) 

6 (t + At) =6 (t) + (e 2 (t) + Q 2 (t)) At + F (t) At 1/2 , (5.8) 



and it is not difficult to show, that the trajectory are continuous everywhere, that is 
9 (t + At) — ► 9 (t), but has no derivative anywhere due to the member ~ At 1 / 2 in (5.8). 

Suppose, that At = t/N, with N — > oo, then eq. (5.7) can be regarded as transition 
probability for 9(t) = 9 k -^ 9 k+ i = 9 (t) at a time At in the model of Brownian motion. 
The eq. (5.7) thus gives the total Fokker-Plank measure of the space Rjqm\. 

Now we can construct the full wave function of random QRHO. Using expression (2.5) 
and turning to moving coordinate system by means of regular shift (4.5) for the wave 
functional the final expression is obtained: 

*+ (n; x, t) = (*+, (n; x, t\9 (t))) = I ' Dy.{0 (t)} *+, (n; x, t\9 (t)) . (5.9) 



{«(*)} 
In (5.9) by Dfi {9 (t)} the measure of functional space Riotm is denoted: 

Dfi {9 (t)} = a- l dfi {9 Q } x dfi {9 t } lim {{2net/Ny- x 

N— >oo 

(5.10) 

N r on 

x n exp -(N/2et){9 k+1 -9 k -(9l + n 2 (t))t/N} 2 d9 k+1 }, 
where a, dfi {9q} and dfi {9t} are determines, accordingly, by the following expressions: 



a= Dfi{9(t)}, (5.11) 

dft {9 } = 8{9 - io (t) /6 (t)) ^o, (5.12) 

d^{0 t } = P(0,t|O,O)d0f (5.13) 

In the formulae (5.9)-(5.13) a is normalization constant for the functional integral 
(5.9) with full Fokker-Plank measure, nonequal to one, integration over dfi {9q} measure 
provides transition to moving coordinate system and integration over dfi {9t} measure 
provides, accordinglythe process of averaging by coordinate distribution 9 in a moment 
of time t. By integration over dfi {$o} measure in expression (5.9) it is possible to get the 
factorization of regular and chaotic motion. Then the wave function of Brownian particle 
will be rewritten as follow: 



V£(n;x,t) = y + (n;x,t)J DjS{0 (*)}*+, (n; a, t|0(i)), (5.14) 

where L>£ {9 (£)} = £>// {# (£)} /d^ {# } • 

6. Solution of the equation for distribution function of stationary 
Markovian process 

Let's consider the probability P(B,t\0, 0) = Q(6,t) that characterizes the distribution of 
the coordinate 9 in the i?{e( t )}-space as a function of time "£". In this case (5.6) should 
be interpreted as the conservation law for probability density 



d t Q (6, t) + do J (9, t) = 0, J (9, t) = - [9 2 + nl(t)) Q (9, t) - edeQ (9, t) , (6.1) 

with the initial and boundary conditions 

lim Q(6,t) = 6(6), lim Q(0,t) = O. (6.2) 

t— ►— oo |6>| — >+oo 

For the boundary fluxes one has Jo = Jo ( — oo, t) = Jo (+oo, t) and it does not vanish 
since (9 2 + Qf)(t)) on that boundaries turns to be infinity. At the limit t — ► +oo the flux 
density turns to its limit value 

J 0/ = lirn^ {j(9,t) sign (9 (t))}, J 0/ = J (^L) • (6-3) 

From equation (4.2) it follows that lim 9 (t) < and as a consequence Jof > 0. 

As a result the equation for probability distribution Q s (9) for stationary process can be 
derived from equations (6.1) and (6.3) 



J 0/ = (9 2 + n 2 out ) Q s + ed e Q s , d e = d/d9. (6.4) 

It may be easily solved, giving 



Q s (e,n out ;9)=e- 1 / 3 Q s (\,r,9)=Jf 5 ex V -- - A 7 J dzexp(j + \>yz) (6.5) 

where A = fan/e 1 / 3 ) 2 , 7 = (^M„) 2 , B = 9/e 1 / 3 . 

The constant Jo/ may be calculated from the normalization condition and has the 
form [13]: 



' .3 



J -/ = ne-^Jo} = tt^V 1 / 3 J dz z' 1 ' 2 exp -^ - X 1Z . (6.6) 

o ^ ' 

For the Jo/ one can obtain another representation via the special functions. It may 
be done by passing to Fourier components in the equation (6.4) [14]: 



J" 1 = ne- 1 / 3 J"/ 1 = 7T5- 1 / 3 



^ 2 (-A 7 ) + Bi 2 (-A 7 ) , (6.7) 



where Ai (x) and Bi (x) are linear independent solutions of Airy equation [15]: 



y — xy = 0. (6.8) 

Numerical calculations of function of distribution Q s (A; 9) = Q s (A, 7 = 1; 9) in de- 
pendence of 9 from (6.5) and (6.6) for some values of parameter A when 7 = 1 are shown 
on fig. 1. 



Fig. 1. Distribution of stationary process Q s ( A; 9) over 9 in dependence of parameter 

1/A~e. 

It is visible, that when e — > 0, i.e. when passing to regular case in initial problem 
(2.1), the function of distribution of stationary process turn to delta- function of Dirac. 



lim Q s I A, 7; 9 

£— >0O,7<0O \ 



' w • 



(6.9) 



7. Calculation of transition amplitude for the random QRHO 

The transition matrix for the random QRHO will be evaluated as a limit t — ► +00 of 
the projection of the total averaged wave function (5.9) on the asymptotic wave function 
^ out (to; x, t) 



t— >+oo 



s mn = l™ (^out(m;x,t)^ r (n;x,t) 



(7-1) 



Taking into account, that measure in the functional integral is real and positively 
defined (5.7) we can change the order of integration in the expression (7.1) and represent 
the transition matrix in the following form: 



qbr 



- lim 



« (*ifc(*)])) {KW]} =£& \ s ™ (t\m) 



mm 



(7.2) 



were S^ (t\[£(t)]) is a stochastic transition matrix, 



S% (t\ m)l) = {*out (m; x, *)*+, (n; s, t| K(t)])) x , [£(*)] = £ (t\W(t)) . (7.3) 

It should be noted, that when Hamiltonian (2.1) is a real function, stochastic matrix 
Smn (*I[£(^)D as wen as ^ s averaged value S^ n are unitary ones. If Hamiltonian is a 
complex function the unitarity of this matrix breaks down. 

Now we will pass to the calculation of expression for transition probability. Taking 
into account completeness and ortogonality of functional basis \I/J C (n;x, t\[£ (t)]) (see 
(3.6)) calculation of stochastic matrix elements S^ (£| [£(£)]) is convenient to carry out by 
generating functionals method. Let us construct generating functional in following form: 

oo n 

*ttc(z,xM(t)}) =E -^Vttc^xMit)]), (7-4) 

n=0 V™ 

where z is some subsidiary complex function. After substitution of expression for 
^ftc ( n 'i x it\[£ (t)]) from (3.5) to (7.4) and carrying out summation [10, 15] we find the 
following equation: 

KAz^Mm=^y^p{-l(ax 2 -2bx + c)y (7.5) 

where a, b and c have following type: 

As it seen from (7.6) the generating functional dependence over the x coordinate is 
stochastic gaussian packet. In a limit t — ► — oo (7.6) turns over to the ordinary gaussian 
packet 

Vf tc (z,x,t\[Z(t)]) -► V in (z,x,t) = 

t — ► — OO 

_ (7-7) 

= (Skjir)* exp {-| (fi in x 2 - 2 v ^~zxe- iQ - t + sV 20 *** + in in t)} . 

The generating function of (out) state can be obtained by making in (7.7) formal 
substitutions Q in — > Q out and z — > z±. 

Now we will consider the following integral: 

I( Zl ,z 2 ;t\[Z(t)}) = (Vtut(zi,x,t)*ttc(z2,x,m(t)})) x - (7-8) 

Substituting expressions (7.5) and (7.9) to (7.8) and carrying out integration by x 
coordinate for the generating functional and generating function of (out) asymptotic space 
one can obtain the following equation: 

/(^,^t|[e(t)]) = (fi m fi oui )^-^yexp|-i^-^|, (7.9) 

where the following notations are made: 

10 



B(t\[£(t)]) = v / 2n~r 1 ^2 + v / 2n^exp(in out t)zi, 



(7.10) 



c (*l [£ (*)]) = ex P (-i2r) z| + exp {i2tt out t) z\ - ifi out t, 

where z = \z\ exp (i arg z) , z — \z\ exp (— % arg 2) . 

As it is easy to ,that the I (zi, z 2 ; t|[£ (t)}) integral is generating functional for matrix 
element S% (t|[£(t)]) ^ ' 



00 m „n 



i(z 1 ,z 2 ; t \iat)]) = i: -±=*=s%(t\e(t))- 

mn=0 VmW. 



(7.11) 



Decomposing / (z±, z 2 ; t|[£ (t)]) into Taylor power series over z\ and z 2 from (7.11) we 
find the following final expression 



21=22=0 



(7.12) 



sS « (*)]) = v^n fa 8 ^ 1 &> *>> *IK (*)])}. 

V m!n! k Jx 

Bellow the expressions for some first stochastic matrix elements are shown without 
some phases irrelevant for scattering process, 



s&(t\[t(t)]) = V2(njn out ) 1 * (-iz/n out + ty l \ sg{t\[z(t)\) = (s£{t\[Z(t)\) 



S02 m (*)])= Si 



stc 
00 



■i + i^V4,fs« ,! 



"OM</ |£|" 



2 V^OO 



exp — i2Vt 



dt' 



(7.13) 



S£(t\[t (!)]) = s- t 



stc 
00 



l + d^Vfc^ 2 



a 



Hill 



^00 



The matrix elements (7.13) after regular shift (4.5) in functional space R{e(t)} have the 
following form: 

S^(t\e(t)) = V2(Q m /Q out y(-t6/n out + iyKxp(-^ J 6(t')dt'Y 



sf\ c (t\e(t)) = (s f (t\e (t))) 3 , 5 2 f (t|e (t)) = 5, 



A'tc 

00 



-1 + 



^nV (cstc\ 2 

nj V 6 °oJ 



, (7.14) 



Sg(t\9(t)) = Si 



stc 
00 



-1 - [g* 



2 /cstc\2 



(5 s * c ) exp -2 / 0(f) dt' 



x 



x exp 



t -2Je(t")dt 

—i2Vt in J e -°° dt' 

—00 



" \ 

) 
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As it visible from expressions (7.13) and (7.14) for stochastic matrix elements the only 
transitions possible are transitions with same evens non-depending from value of fluctu- 
ation constant e. As to symmetry of matrix elements by oscillation quantum numbers of 
initial and final channels n and m, it breaks down as a result of irreversible character of 
quantum mechanics constructed here. 

To demonstrate the proposed approach we shall represent the results of evaluation 
of the " vacuum- vacuum" transition probability under the condition of Wiener's process. 
Using equations (5.9)-(5.11) and (6.5)-(6.6) from (7.1) one obtains: 



S* (A, p) = (1 - p)Hh (A, 7 ) - ih (A, 7 )} , 



(7.15) 



MA.")= / deU^Q s (\, T ; 



(7.16) 



j 2 (A, 7 ) = / de-J d —^Q s {\r,e), rf(A, 7; ^ = [i + — 



7T2 \ 2 



(7.17) 



Here p is a reflection coefficient of the correspondent one-dimensional quantum prob- 
lem (see [16]), 7 (p) is denoted by the barrier shape i.e. by the frequency fio(t). I n the 
case of the step-shape (fig. 2) barrier one has: 



l{p) 



O ^ 2 

J L out 



1+p 1 ^ 
1- P V2 



(7.18) 



Fig. 2. Barrier model of dependence of frequency Q over time t. It is clear, that changes 
of S in range of < 5 < oo cause changes of reflection coefficient p in range of < p < 1. 

As a result, using equations (6.5)-(6.7) and (7.14)-(7.17), for the probability of "vacuum- 
vacuum" transition one obtains: 



Ato(A,p)= S^(X,p) 



(7.19) 
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^(A,p)| 2 = V / W{/i(A,7) + /|(A,7)}. 



(7.20) 



The result of calculation of the transition probability (7.18)-(7.19) are represented on 
(fig. 3-4) as a function of p and A. 



Fig. 3. " Vacuum- vacuum" transition probability in dependence of A and p. 



Fig. 4. Dependence of "vacuum- vacuum" transition probability over A in the case when 

p = 0. 

As it visible from fig. 3. the probability of " vacuum- vacuum" transition starting from 
some value of A (or e) have nonmonotonic behavior in depending over reflection coefficient 
p. This fact distinguish the stochastic problem from regular one. 
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8. The vacuum thermodynamics for the asymptotic space S s cH 

It is well known, that the principal object of interest for quantum statistical mechanics 
is the density matrix p(x,x') , that after Dirac and von Neyman (see, [17]) is determined 
by expression 



p{x,x') =Y^ Pk^k(x)(p k (x') : (8.1) 

k 

with distribution function for canonical distribution P^ = exp (—/3Ek), and wave func- 
tion being the solution of Schrodinger equation Hip^ = E^ipk, f3 = {kT)~ with T being 
the temperature of the system and k being the Boltzmann constant. 

Since for the our problem the wave function of the system is the complex stochastic 
process it is natural to turn to thermodynamic description. It is possible to develope such 
picture even for the group of states, corresponding to quantum number "n". 

Here we discuss the thermodynamics of vacuum in the asymptotic space 

E s = R l as <g> R {Ws} = t Hm R 1 <g> R{ W (t)}- 
Definition l.The stochastic density matrix for vacuum in the space S s 

p stc (x, t; 9{t) | x', *'; 6(f)) = {*+ c (0; x, t | #(£)) ^+ c (0; x',£' | 0(t'))} . (8.2) 

Definition 2. The expectation value of stochastic operator A(x,t \ 0(t)) in the vacuum 
have a next form 

( A \ ac =JSSo { Tr * (<^-*> w J /Tr * (<**W } ' ( 8 - 3 ) 

with being the trace over the coordinate x. 

Definition 3. The nonequlibrium partition function of the vacuum and quantum oscil- 
lator system is 

§Z c {e,n as] t) =Tr x {(p stc ) {m} } . (8.4) 

Knowing the partition Function it is easy to determine all the thermodynamic prop- 
erties of the system: 

a) average internal energy 



U mc (e, tt as ) = Jim U mc (e, tt as ; t) , U„ ac (e, fi as ) = -d e {lntf^ (e, fi as ; £)} , (8.5) 

t — >+oo 

b) free Helmgoltz energy 

F vac (e, tt as ) = -e' 1 Jim {\wd vac (e,Q as ;t)} , (8.6) 



c) the Entropy 

The practical computations start from the stochastic density matrix 
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p stc (x,f,e(t) | x',t'-6(t')) = (n as /n)i exp{-n as (x 2 + x ,2 )/2- 

(8.8) 
-\ J 0{T)dr-\ ) e(r)dr-i[e(t)x 2 -e(t')x' 2 ]}, 



with Q as being the frequency in the asymptotical space S s .The stochastic Hamiltonian 
(2.1) is discussed as an example of stochastic operator A(x,t \ 0{t)). After the nondif- 
ficult calculation from (8.3) with taking into account (8.8) one can obtain the energy of 
vacuum+oscillator system 

oo 

E (A; n as ) = -±n as J of J dzz-^ 2 exp (-^ 3 /12 - Xz) + 

o 

+ \n as |l - ij 0/ 7 dzz^ 2 exp {-z*/l2 - Xz) + (F) /4^\ + ( 8 . 9 ) 

oo 

+i^tt as J of J dzz 1 ' 2 exp (-^ 3 /12 - Xz) . 

As it clear from (8.9), the first term in the energy expression diverges, corresponding to 
the infinite energy of the vacuum. The second term corresponds to the oscillator energy, 
that is shifted by the interaction with vacuum 

E™ (A; n as ) = hl as {l~\d a [(1 + a) d a \nA(-X + a)] | Q=0 } , (8.10) 

A (-A + a) = At 2 (-A + a) + Bi 2 (-A + a) , 

and <F> in (8.10) in our situation is zero. Note, that second term in (8.10) is 
analog the Lamb shift of the energy level it is well-known from the standard quantum 
electrodynamics [18]. The third term in (8.9) corresponds to the width of the energy 
ground state and is inverse proportional to its decay time 

a/A 
At = 2— {d a \nA(-X + a)}\ a=0 . (8.11) 

Now let us calculate the basic thermodynamical function of the system - the Entropy 

We start from the expression for partition function (8.4) 

$ vac (e;t) =Tr x {( Pstc ) m } = BoWBfatyt), (8.12) 

with B (t) and Bi(e,Q;t) being determined by the conditions 

£o(t) = (exp (- f e(t')dt'\\ = f d9u(9,t), (8.13) 



{«(*)} 



— oo 



B l (e, Vt as ; 9,t)= J d9Q (e, Vt as ; 6, t) , (8.14) 



15 



and u (9, t) in (8.14) by Feynman-Kac theorem being the solution of parabolic equation 

[19] 

d t u(9,t) = ^d 2 e u(9,t)-9u(9,t), (8.15) 

with initially and boundary conditions of the type (6.2). 

It is easy to show, that those solutions do not depend upon the volume of e and the 
limit of Bo(t) at t — > +oo is 2 -1 / 3 . But the function B\ (e,t) is clearly dependent on e. 
The general solution of (6.1) with additional conditions (6.2) can be represented in the 
form 



Q(e,n as ;9,t-t')=Y / e- x ^ t -^Q k s (e,n as ;9), t y t' , t' -► -oo, (8.16) 

fc=0 

with A = 0, and Q° s (e, fi as ] 0) = Qs{\ 1\&). So after differentiating of (8.12) and 
taking into account (8.14)-(8.16) in the limit of t — ► +oo one gets for the average internal 
energy the expression 

oo 

U mc (e, Q as ) = - f d9d £ Q s (e, Q as ; 9) . (8.17) 



After straightforward computation we have 

U mc (e, n as ) = U mc (A) = ±- {1 + 2\d a (In A(-X + a)) \ a=0 } . (8.18) 

Taking into account (8.12)-(8.15) it is possible to have the expression also for Helmholtz 
Free Energy 

F vac (e,Q as ) = ^-\n2, (8.19) 

and for the Entropy of the vacuum can be represented in the form 

S vac (e, Q as ) = S vac (A) = — {8 a In A(-X + a)} \ a=0 +- (1 - In 2) . (8.20) 

That give the expressions for all thermodynamical potentials of quantum oscillator in 
the ground state, interacting with vacuum. The fig. 5. shows dependence of energy of 
oscillator "ground state", its shift and entropy of vacuum over parameter A in units of 
Boltzmann constant k. It is visible, that when system turn to balance state (e — > 0, i.e. 
A — ► oo), the entropy aspire to maximum value. 
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9. Conclusion 

Chaos in quantum systems was observed first by one of the founders of quantum mechan- 
ics Wigner, when he studied the nucleus energy spectrum [4]. However, Wigner and many 
other researchers associated this phenomena with nonclear and exotic nature of nucleus 
interactions, and they have exit from this difficult situation by introducing some urge 
parameters to the nucleus statistical theory. But, as it was shown by further investiga- 
tions, the chaos arise in spectrum of quantum systems with some particular interaction 
potentials, for example, hydrogen atom in strong magnetic field. Some modern researches 
of modeling of bimolecular chemical reaction [5-6] showed, that chaos affects the wave 
function of quantum system. In other words, we have obvious example of violation of 
deterministic principle related to the basic object of quantum mechanics - wave func- 
tion. To overcome this difficulty the authors in the framework of internal time idea (see 
[9], [20]) have developed new representation for multichannel scattering. It was proved, 
that system, including three or more particles, in general case have chaotic internal time. 
The last shows, that constructed quantum theory in general case being irreversible in 
relation to that time. It must be noticed, that chaos may be caused not only by the 
difficult dynamics of the quantum system, but also by the strong interaction of system 
with thermostat (with vacuum in our case). This situation was investigated in framework 
of one-dimensional random QHRO model. The main idea consist in representation of the 
wave function as a complex probabilistic process ty s t c {x,t\W(t)) on the extended space 
H = R l <g> Rswit)}- Using the model one-dimensional nonlinear Langevin equation the 
separation of variables in initial SDE for wave function was made and stochastic basis set 
tyf tc (n; x,t\W(t)) of quantum system was obtained. One of the very important features of 
such representation is that at least for the closed system "vacuum+oscillator" the nonlin- 
ear Langevin SDE generate real full Fokker-Plank measure in the functional space Rswit)}- 
This circumstance provide exact mathematical basis of the constructed mixed functional- 
wave representation of random QHRO wave function \l/ br (n;x,t). The developed theory 
unificate two inconsistent concepts: the quantum analog of Arnolds transformation, that 
don't admit arising of chaos inside the trajectories beam,that described by the similar 
topology, and functional integral method, that allows to run over the current tubes of 
arbitrary topology in functional space Rswit)} an d to generate chaos. In another words, 
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the proposed theory allows to establish the connection between the chaotic classical and 
chaotic quantum regions. In this work for the case of Wiener measure exact expression 
for the amplitude of " vacuum- vacuum" transition probability A^ (X, p) was constructed 
and it was shown, that behavior of this probability by changing the reflection coefficient 
p of one-dimensional quantum problem is nonmonotonic. It is in detail investigated the 
properties of " vacuum+oscillator" system in asymptotic space E s = R^^R^Ws} an d it was 
shown, that the ground state of oscillator is described by innumerable basis set in Hilbert 
space, unlike the case of standard quantum mechanics. The thermodynamics of vacuum 
in asymptotic space S s is studied in detail and the energy of oscillators ground state with 
analog of Lamb shift and level width calculated, the expression for the Entropy of system 
in dependence of couple constant A is constructed. Let us remind, that the Lamb shift 
of energy levels in hydrogen atom from the point of view of quantum electrodynamics 
is obtained in framework of perturbation theory, but in proposed theory the analog of 
Lamb shift is obtained without including of perturbation theory. The last feature of this 
approach indicate its nonperturbative nature, as was expected early [11]. And in the 
end let pay attention to the principle difference of this theory from any other quantum 
approach, that is the possibility of decay of ground state. 
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